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magnetic field were measured at a land site (La Mesa) and 
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I. INTRODUCTION 


This thesis research is part of an ongoing effort being 
conducted by the Naval Postgraduate School to obtain a long- 
term data base of ELF electromagnetic noise on the ocean 
floor. The result of this project should be a better under- 
Standing and interpretation of that noise and its sources. 

Magnetic noise iS important to the geophysicist and to 
the Navy. Magnetic techniques are applied in exploration 
for minerals, the prediction of communications propagation 
and to the modeling of the sphere on which we live. Geo- 
magnetic studies are of interest to the Navy for their 
application to submarine detection, submarine communications 
and to mine warfare. 

The objective of this thesis within the broader context 
of the overall project was to design and validate a computer 
code that would facilitate data manipulation and analysis. 
To that end, the present computer code calculates for two 
Sites (one on the ocean floor, the other on land), the 
individual site's power spectral densities and Stokes 
Meweametcers 0, 1, 2 and 3 for two of the individual site's 
orthogonal components in the horizontal plane. The code 
also computes the coherence of these two orthogonal compo- 
nents and their phase difference. Lastly, the program 


computes the right and left circular polarization power 
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Spectral densities of the sites and the coherence between 
the seafloor and land-based sites. 
As discussed, the present code applies only to two 
Sites and two orthogonal geomagnetic components at each 
site. However, extension to several sites and/or components 
Should be an easy undertaking due to the style of the present 


program, 
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II. BACKGROUND 


A. THE MAIN GEOMAGNETIC FIELD 

Historical observation and modern measurement demonstrate 
that the principal source of the geomagnetic field (approxi- 
mately 50,000 nanoteslas) is beneath the surface of the 
earth. An easy-to-picture causal interpretation of this 
phenomenon is painted by imagining the earth with a dipole 
magnet near its center as pictured in Figure 2.1. The 
magnet axis 1S tilted at an angle of 11.5 degrees from the 
earth's rotation axis. This model provides a prediction 
accuracy to a few percent of the measured geomagnetic field. 
Figures accurate to less than one percent of the earth's 
measured field can be obtained uSing more sophisticated 
models [Ref. 1]. 

The generally accepted physical process for the genera- 
tion of the geomagnetic field is that of the self-sustaining 
dynamo [Ref. 2]. Figure 2.2 presents a simple picture of 
such a dynamo. The physical requirements of this theory 
are mostly met in the interior of the earth. For example, 
the earth's center is a large conducting system. Addition- 
ally, we can assume that the earth's molten core is of a 
differential rotational nature with turbulent upwelling. 

Energy sources for such a dynamo are most likely avail- 


able in the earth's core due to a number of processes. This 


tes 








Eccentric Dipole Model of Geomagnetic Field 


Prugure 2. 1 
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Axis of rotation 







Nonrotating coil 


Rotating disc 


Figure 2.2 Simple Disk Dynamo 


all boils down to mean that the motion of a conductor, such 
as the earth's molten core, in a magnetic field produces a 
current which in turn produces a magnetic field in support 
of the original magnetic field [Ref. 3]. 

A Simplified picture of the dynamo theory originates 
by assuming that in the distant past a magnetic "seed" field 
was captured by the earth's conductive core. The differen- 
treroaation Of this core "winds up" the field into an 
intense azimuthal field, such as would happen to a twisted 


rubber band. Turbulent upwelling of the core tends to 


ie, 





carry this azimuthal field outward, converting it to a 
toroidal helical field within an associated azimuthal 
emcimating Current. This current in turn generates the 
earth's dipolar field. 

An additional advantage of the self-sustaining dynamo 
theory lies in its ability to explain the well-documented 
reversals of the geomagnetic field. It has been postulated 
that two or more coupled "dynamos" as shown in Figure 2.3 


can cause reversals of the main geomagnetic field over 


geologic time scales [Ref. 4]. 





Pigure 2.3 Twin Disk Dynamo 


ES 





B. GEOMAGNETIC FLUCTUATIONS 

The geomagnetic field changes with time. The changes in 
the field have periods ranging from fractions of seconds to 
thousands of years. The very slow variations (periods of 
years Or more) are not of interest here but are known as 
secular variations and are thought to be geologic in origin. 

The geomagnetic field changes of interest here are those 
fluctuations that occur over the period of a day or less. 

It 1S customary to view the magnetic field as the stable 
main field with these comparatively rapid fluctuations 
Superimposed on it. Figure 2.4 presents the spectrum of 
the superimposed fluctuations. 

The fluctuations (Figure 2.4) that are superimposed on 
the main field range in intensity from .05 to 50 nanoteslas 
and have various causes. Some major contributors to these 
fluctuations are in the ionospheric currents overhead and 
in magnetic storms which probably originate in solar 
eruptions. 

Finally, arriving at the point of interest here, are 
those disturbances to the earth's magnetic field that 
propagate through the earth's magnetosphere as hydrodynamic 
waves. These small perturbations cover a wide range in 
amplitude and duration and are known as micropulsations. 
Figure 2.5 presents a graphical display of these disturbances. 


Micropulsations are thought to be the result of fluctuations 
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in the incident solar wind affecting wave-particle interaction 
in the magnetosphere [Ref. 5]. As they propagate through 
the conductive layers of the atmosphere, they are trans- 
formed into electromagnetic waves via ionospheric currents, 
which we detect at the earth's surface as magnetic field 
fluctuations. There are two major categories of micro- 
pulsations: continuous and irregular. 

Smoothly varying, periodic pulsations are classified as 
(Pc) micropulsations. Their distinct periodicity is indica- 
tive of a resonance inside the magnetospheric cavity. One 


model proposed to explain this occurrence of smooth, periodic 
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Pulses 1S that ‘the resonance is a result of an oscillatory 
hydromagnetic wave. This wave is associated with groups 
of ions which "bounce" back and forth in the magnetosphere 
[Ref. 6]. 

iiese COntinuous™micropulsations are further classified 
according to their periods as follows: 


ieeect:. (0.15-5%second period) 
Known as "Pearls," these micropulsations 
are generated by the cyclotron instability 
of the energetic protons. They have been posi- 
tively correlated with solar disturbances and 
Occur during daylight hours in the auroral 
zone and during night and early morning hours in 
the midlatitudes. Their average amplitude is 
approximately 0.05-0.1 nanotesia [Ref. 7]. 


pee c2: (5-10 second period) 
This is a diurnal phenomenon that shows some 
positive correlation with solar activity and the 
seasons. They usually decrease in their period 
as Magnetic activity increases. Their average 
slteude 1590. 1-1 nanotesia. 


3. Pc3: (10-45 second period) 
This micropulsation is similar to Pc2 pulsations 
except for the difference in period. 


4. Pc4: (45-150 second period) 
Sunspot activity appears to have an effect on 
Pc4 pulsations. The Pc4 pulsations also change 
their period seasonally with an average amplitude 
of 5-10 nanotesla. 


Be ©’ G5: (150-600 second period) 
These large-scale pulsations typically occur 
during morning and evening with amplitudes of 
10-100 nanotesla. Their duration shows a strong 
geomagnetic latitude dependence. 
Irregular (Pi) micropulsations appear to be associated 
With polar magnetic disturbances and are observed during 


night-time hours. The theory is that these micropulsations 


Ze 





arise from disturbances of the ionosphere due to particle 
bombardment from the magnetosph2re. Increased solar radia- 
tion levels also appear to have an effect on their existence. 
Irregular micropulsations are classified as follows: 
1. Pil: (1-40 second period) 
These pulsations have a small amplitude and 
usually occur at night and early morning. 
They vary in intensity from 0.01-0.1 nanotesla. 
They demonstrate a positive correlation with 
auroral disturbances. 
Ze Pi2: (40-150 second period) 
Pi2 pulsation amplitude range is on the order 
of 1-5 nanotesla. They usually occur during 
the early evening hours but may continue through- 
out the night. The periods of these pulsations 
decreases with increasing magnetic activity. 

The Pc5 and Pc4 micropulsations last several hours. 
However, the Pc5-1 micropulsations and irregular pulsations 
have a maximum duration of approximately one hour but may 
cease after only a few minutes [Ref. 8]. Above 3 Hz, the 


most pronounced field variations are the earth ionosphere 


cavity resonances of Schumann resonances [Ref. 9]. 


C. OCEAN WAVE EFFECTS 

In addition to all of the geomagnetic fluctuations 
discussed above, if working in the ocean environment, one 
must take into account hydrodynamically induced magnetic 
fluctuations. These fluctuations are internaily produced 
resulting from the movement of a conductor (seawater) in 
an external magnetic field. The effects of waves have been 


calculated [Ref. 10] and are pictured in Figure 2.6. 
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Weaver's equations have been modified to incorporate bottom 
reflection and transmission with the result that the effect 
of bottom decreased the magnitude of the induced magnetic 

field by less than 3 dB for frequencies greater than 0.16 Hz 


in water depths less than 500 meters [Ref. 11]. 
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Ill. COLLECTION SYSTEM 


A. DATA COLLECTION SYSTEM 
The data acquisition system is illustrated in Figure 3.1. 

Major components consist of: 

(1) coil antenna sensors (2) 

(2) preamplifiers (2) 

(3) signal conditioners (2) 

(4) pulse code modulation system (1) 

if) radio transmitter (1) 

(6) radio receiver (1) 

(7) tape recorder (1) 
For the work reported in this paper, two systems consisting 
of components (1) through (5) were deployed. They were 
Similar in most respects with the major exception that the 
ocean equipment was enclosed in watertight glass spheres. 
A thorough description of the data acquisition system is 
given in Reference 12. 

1. Coil Antenna Sensors 
Each coil antenna sensor consists of two continu- 

ously wound coils which are constructed of about 5460 turns 
of 18-gauge copper magnet wire. Glass spheres are used to 
maintain the watertight integrity of the coils and all other 
immersed equipment for the ocean configuration. See Figure 


3.2 for the ocean system. The coils are mounted orthogonally 
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in a nonmagnetic frame with each coil having a resistance 
of 120 ohms and an inductance of approximately 9.31 Henries. 
2. Preamplifiers 
The preamplifiers are model 13-10A low-noise ELF 
amplifiers manufactured by Dr. Alan Philips. The final 
Stage of the amplifier contains an active low-pass filter 
which provides a sharp cutoff for frequencies above 20 Hz. 
Frequencies above 10 Hz are not examined in this paper. The 
overall amplifier gain for inputs of less than 2.5 milli- 
volts is 60 dB in power. The purpose of the filter is to 
remove 60 Hz noise and to prefilter the data before 
eet zation. 
See Olgnal Conditioners 
The signal conditioners receive analog signals from 
the preamplifiers, amplify them by about 30 dB and limit 
Signals with peak amplitudes of 7.5 volts (modifiable) from 
entering the pulse code modulation system. 
4. Pulse Code Modulation System 
The pulse code modulation system was designed and 
manufactured by Dr. Robert Lowe, Lowecom, Inc. The system 
offers a 15-channel analog input capability with selectable 
Sampling rates of 2, 4, 8, 16, 32, 64 and 128 samples per 
second. The PCM system uses a cryStal oscillator and inte- 
grated circuitry to develop the clocking pulses, a 16-channel 
analog multiplexer, a 16-channel, 12-bit analog-to-digital 


converter as the basis for the pulse coding. The crystal 


ao 


| 3 








clock oscillator loses 1 bit in 1 million. The clock pulse 
gates the analog multiplexer, the analog-to-digital converter 
and additional circuitry that form the pulse code words. The 
PCM output is organized in frames with each frame headed by 
a synchronization word which is followed by the pulse-coded 
samples from PCM channels 1 through 15. 

5. Transmission and Recording 

After encoding, the data is transmitted via VHF 

radio to a Naval Postgraduate School located receiver and 
recorded on an instrumentation tape recorder. For the data 


of this paper, a sample rate of 32 samples per second was used. 


B. DATA ANALYSIS EQUIPMENT 
1. Analysis Hardware 
The PCM data is fed into a PCM decoder which places 
the decoded PCM onto a 9-track 800 BPI computer tape. This 
computer tape 1S read into the Naval Postgraduate School IBM 
3033 computer for manipulation and analysis. 
2. Data Analysis Software 
The IBM 3033 is used to convert the data to usable 
form and then to perform any operations on the data as desired. 
At present the computer code is written in Fortran IV and 
converts the digital data to voltage, Fast Fourier trans- 
forms the data, calculates power spectral densities, deter- 
mines Stokes parameters 0 through 3 and provides phase and 


coherence figures. Appendix A contains test verification 
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of the code. Appendix B contains the code that prepares 
data from one site for storage in the IBM 3033 Mass Storage 
System for later recall and comparison with another site. 
Appendix C contains the code that prepares and manipulates 
the code from a second site, retrieves from the Mass Storage 
System that data previously prepared and manipulates the 


two data sets simultaneously. 
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Vie COMPUTER SCODE 


A. INTRODUCTION 

As mentioned earlier, one of the intentions of this 
effort was to develop a code that would permit the comparison 
and/or simultaneous geomagnetic analysis of two different 
Sites. The present data collection system provides for the 
recording of transmitted data from those sites onto different 
tracks of the instrumentation tape recorder. Since only one 
track of that analog tape can be placed on one digital 
computer tape, a means of storing one channel had to be 
devised. As can be seen by a quick review of Appendices B 
and C, the present technique that allows this simultaneous 
manipulation of data from two sites is to read the computer 
tape from one site and to process that data to the point 
where it 1S ready for analysis. These results are then stored 
in the IBM 3033. Subsequently, the same routine (with modifi- 
cations to fit the parameters of the second site) are applied 
to data read from the second site computer tape. Now both 


sets of data are in the same form and ready for manipulation. 


B. DATA ANALYSIS SOFTWARE 

The two programs (one for each of two sites) are 
essentially the same with the exception that one program, 
Beooream 6B. isea short version of the other, program A. Thus, 


the next paragraphs will discuss in detail the longer program. 
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This description fits equally well to program B except that 


program B ends after applying the system transfer functions 


and placing these results in the Mass Storage System. 


1. The Main Program 


The main program 1s divided into sections which 


perform the following functions: 


*& 


* 


Data input 

Fourier analysis of time series data 

Application of system transfer function. Note: At 
this point, program B reads the data into the Mass 
Storage System and ends. Program A, at this point, 
retrieves that same data from Mass Storage. 

Data averaging 


Gateulation and plotting of: 


(1) power spectral densities of the individual 
Site geomagnetic orthogonal components. 


(2) orthogonal component coherence of the geo- 
magnetic field of one site. 


(3) Stokes parameters 0 through 3 of the individual 
Site orthogonal components. 


femoinele site right and left circular polarization 
power spectral densities. 


(5) single plane coherence of right and left 
circular polarization between sites. 


a. Data Input 


The data is read into the IBM 3033 via the sub- 


routine "read'' which was written by Tim Stanton of the 


Naval Postgraduate School. The "read" subroutine reads tue 


PCM 


data from the computer tape and converts it into an 
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imeeecer trom 0 to 4095, The main program than converts 
this integer to voltage. See Appendix B or C., 
b. Fourier Analysis 
The next section of the code utilizes a sub- 
routine developed at the Naval Postgraduate School and 
oe d “rourt.’  "“Fourt" Fast Fourier transforms the voltage 
data from the time domain to the frequency domain. See 
Appendix A for validation of this routine. 
eee transfer Functions 
After the Fast Fourier transform, the program 
applies the transfer functions for the system. Development 
of the transfer functions is contained in Reference 12. 
The transfer functions are applied up to 10 Hz as straight 
line functions. The data is now in nanoteslas. 
d. Power Spectral Density 
This quantity is computed by multiplying the 
previously calculated spectral density summation by the 
observation period and dividing by the number of data 
groups (typically 20) that the computer tape was segmented 
into, providing an average power. 
eee Cone rence 
The next steps of the program determined for one 
Site the geomagnetic orthogonal components Stokes parameters 
0, 1, 2, and 3, and the phase of the components and their 


coherence. Next, right and left circular polarization 





parameters are computed for the individual sites followed 
by calculation of the coherence of these parameters between 


the two sites. 


GeeeeaRAMETER THEORY 
1. Stokes Parameters 
A useful tool for determining the state of polariza- 
tion of a wave are the Stokes parameters 0, 1, 2, and 3 
[Ref. 13]. 
Two linearly polarized plane waves with field 


and €. can be described as: 


vectors in the directions of Ey 7 


# .2 gp , (i Kex - iut) 


eps) hea toe) 


These waves can be combined to give a general plane wave 


propagating in the direction k = kn: 


n~ 


B(k,t) = (2, E, + @, E,) ee 1) 


ea Ey and E, 


polarized wave with its polarization vector making an angle 


have the same phase, (1) represents a linearly 


9 = tan”? (E,/E,) with E, and with a magnitude E = vE, + ES, 
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Use Ey and E, have different phases, the wave (1) is 
elliptically polarized. 

In the simplest case, circular polarization, Ey and 
E, have the same amplitude but differ in phase by 90°. The 


Wave (1) is then: 


E (x,t) = E, (€, * &,) e (k-x - iwt) (2) 
Then the components of (2) are: 
-> 
EY (x,t) = ES COS. (kz = ut) 
() 


Ey (x,t) = E, sin (kz - wt) 


Thus, at a fixed point in space, the fields (3) are swept 
around circularly at the frequency w. For (E, ie the 
rotation is counterclockwise if the observer is facing the 
oncoming wave. This wave is described as having positive 
helicity. (Ey - eS.) has negative helicity or is rotating 
clockwise when looking into the wave. 
The two circularly polarized waves (2) form an 

equally acceptable set of basis fields for description of 
a general state of polarization by introducing the complex 


orthogonal unit vectors: 











ae P% a 
€,= ‘See ge 
Pet? 
Suen that 
ax n = 
oan ° ae = 0 
x a 
5° 0 
AX an 
oa ° oan = 1] x, 


Then a general representation of (1) is: 


>> 
EGE (Ee Ee) eg ee 2 WE) a 
where E, and E| are complex amplitudes. If E, and E_ have 
different amplitudes but the same phase, (4) is an elliptically 
polarized wave with principal axes of the ellipse in the plane 
defined by =F and eae 

Thus, the polarization content of a plane wave is 
known if it can be written in the form (1) or (4) with known 
coefficients (E,, E,) one (GE E_). Conversely, how can we 
determine from observations on a wave the state of its 
polarization? The answer lies in the use of the four Stokes 
parameters. Their measurement or calculation determines 
the state of polarization of a wave. 


For a wave propagating in the z directicn, the scalar 


PGOauc ts : 


i 


<> _ 
er — 
iS a= 





are the amplitudes of radiation respectively, with linear 
polarization in the x direction, linear polarization in the 
y direction, positive and negative helicity. The squares of 
these amplitudes provide a measure of the intensity of each 
type of polarization. Phase information is obtained from 
cross products. If we assume the scalar coefficients of 


(1) and (4) to be a magnitude times a phase factor: 


16+ 16- 
e e 


and use the circular polarization basis (€,, €_), we have: 


s = |, - Bl? + |é. + El? = a2 + a? 
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peeeeeome rence 
If a plane wave of frequency v is propagating in the 


geotrection, it has components: 


Em Ge) a, (t), 79! - L2mvt 


B, (t) S ape) * Pah = i 20vt 


that are mutually orthogonal and at right angles to the 


direction of propagation [Ref. 13]. Then the complex 


degree of coherence can be written as: 


2X yy 
where Jortee are the elements of the coherency matrix 
. : 191-162 
EEL. EE. | a4 a,a,e 
a : 
- = 1$,-1id¢2 
ae EYE, a,a,e a. 


The complex factor Hoy is the complex degree cf coherence and 
is a measure of the correlation between components of the 


field vector in the x and y directions. 
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D. SOFTWARE TESTING 

Miembasis Set for the code calculations is that of a 
circularly polarized wave. Thus a circularly polarized 
wave was chosen as the test input. The wave was created by 
generating a 4 Hz sine wave for sites A and B x coil input 
and 4 Hz cosine wave for sites A and B y coil input. The 
results of this test are valid only at the input frequency 
as a measure of the code veracity. One of the characteris- 
tics of any FFT routine is that small residuals are created 
throughout the spectrum of interest. These residuals are 
primarily caused by the finite observation time. These 
small-valued artifacts will be carried throughout all subse- 
quent parameter calculations and will appear in the test 
peers. This inclusion of unexpected values at frequencies 
Other than the test input of 4 Hz should be considered as 
the result of the artificiality of the single value test 
input. Input of real world data that covers a full spectrum 
of values obscures the "FOURT" routine propensity to create 
these artifacts. 


1. Sites A and B Orthogonal Component 
Power Spectral Density Test 


Spectral density is defined as: 
eilingss * 
“pam 1 1 Cra (05 D 


for site A x coil where x is the Fourier Transform of some 
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input to the x coil. It follows that the average power 
emectral density is: 
Time 


n #e 
PSD 4 = x ; (X,) (X) = a - Time 





where time is the observation time. Thus, for the test 
inputs of sine 4 Hz into the x coils at sites A and B and 
the cosine 4 Hz into the y coils, we expect and obtain (see 
Figures A.1, A.2, A.11, A.12) a peak magnitude of .5 at 
4 Hz at each site. 

2. Sites A and B Orthogonal Component Coherence Test 


Coherence as defined above is: 


for two orthogonal components where eso is defined by: 


. n 
_ Time * 
oe N , (X,) CS) 





S 
xya 


With X and Y our transformed inputs and considering the same 
inputs, we expect to see a coherence of magnitude 1 at a 
frequency of 4 Hz and its harmonics. Figures A.8 and A.18 


demonstrate the validity of the coherence calculations of 
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the code. Also the value of 1 for the test case coherence 
at frequencies other than 4 Hz is a manifestation of the 
residual effect discussed above. 

3. Sites A and B Stokes Parameter 1 Test 


Using the definitions developed in preceding 


paragraphs, Stokes 1 is: 


I 
>) 


Sia pecomerscinets. = 6,) = 2[S, - S| 
for = = Sy as is the case for this test input and as is 
shown in Figures A.4 and A.13. 
4. Sites A and B Stokes Parameter 2 Test 
Again using definitions developed in the preceding 


paragraphs, Stokes parameter 2 is: 
5 = 2a a sin Cee ort UP sonea le" = 6) = 0 


2a 


Figures A.5 and A.15 validate the code for this calculation. 
S. Sites A and B Stokes Parameter 3 Test 
Stokes 3 1S: 


Speeceas =a> =) -(.707)° -(3207)° 2-1 


Figures A.6 and A.16 validate the code for this calculation. 
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6. Sites A and B Stokes Parameter 0 Test 


Stokes 0 is: 


(See Figures A.3 and A.13.) 
7. Sites A and B Phase Difference of x and y Test 
The phase difference of sine 4 Hz versus cosine 4 Hz 
is 90° or 1.57 radians. (See Figures A.7 and A.17.) 
8. Sites A and B Two Orthogonal Component Power 
Spectral Density with Right and Left Circular 
Oolarization Basis Set Test 
The development of the spectral densities for this 
case is similar to the development of spectral density for 


a Single component. Starting with the cross spectra density 


definition for right and left circular polarization 


os, 
a2) 
I) 


(X + 1Y) (right polarized) 


bo 
= 
7 


(X - 1Y) (left polarized) 


Mamrollews that the spectral density is: 


a es x 
aoa ay : (BP ya) (BP ya) (right polarized) 
1S (BM BM) left polarized) 
—— Z 
SM ya q : (B aye xyae (left polarize 
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The power spectral densities will be: 


. Th 
sp = Lime 5 sp 
xya N 1 xya 
. nN 
sM ss = Tame 5 sm 





where the time is the observation time. Examples of this 
calculation with a test input of 4 Hz right circularly 
polarized wave are in Figures A.9, A.10, A.19 and A.20. 


9. Coherence of Sites A and B Right and Left Circular 
Polarization Basis Set Test 


Right and left circular polarizations are: 


SP 
COP = BES (right polarization) 
ve ob 


COM = (left polarization) 


where SP. and SM, are the cross spectra densities of 


bxy bxy 
the two orthogonal components at the sites A and B. 


: n 
_ Time 
SP = : ae 


abxy N (right polarization) 
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_ Time 3 , 
ea xy 7 ; CEM ya? My? (left polarization) 





The test case of sine 4 Hz as the xa and xb inputs versus 
cosine 4 Hz as the ya and yb inputs yields a coherence of 0 
for left circular polarization and 1 for right circular 


polarization. (See Figures A.20 and A.21.) 
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V. ANALYSIS 


Preliminary, "eyeball" analysis was conducted on data 
Semerectcd OVer a period from 1122 local 17 Aug 82 to 0639 
local 18 Aug 82. The timeframes inspected within this 


Window were: 





TIME (local) DATE 

1122-1254 17 Aug 82 
1301-1408 17 Aug 82 
2028-2200 17 As 82 
0121-0251 18 Aug 82 
0259-0431 18 Aug 82 
0507-0639 18 Aug 82 


Original intentions were to collect data over a full 24-hour 
period. However, the batteries at the Monterey Bay site went 
dead after 20 hours. Of the 20 hours of collected data, 
elght hours appeared dominated by system noise and were not 
included in the analysis. However, the remaining data of 
approximately eight hours spans a 24-hour period and provides 
evidence of some obvious trends that lend themselves to the 
cursory analysis conducted here. Representative plots of 

the calculated parameters are in Appendix D. Technical 
Report NPS 61-83-005 contains a complete set of plots of the 
calculations that span the full 20 hours of data coilection. 


This complete set of plots was used in the analysis that follows. 
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ONE SPECTRAL DENSITIES 

Both sites exhibit a marked decrease in power spectral 
density with frequency. There is an obvious exception to 
this observation at 4 Hz at the land site. The sharp peak 
at 4 Hz is attributed to 60 Hz aliasing which is a result 
of the Fast Fourier Transform (FFT) subroutine. Fast 
Fourier algorithm construction is such that data being 
transformed may be convolved back onto itself [Ref. 15]. 
This convolution 1S routinely avoided by using a sample 
rate that is twice the upper frequency of the domain of 
interest. A sample rate of 32 samples per second was used 
in this experiment and should have precluded aliasing of 
the frequency of 60 Hz, especially when combined with the 
amplifier filter cutoff of 20 Hz. However, the 60 Hz 
component of the magnetic field was measured and found to 
be 4 nanotesla/Hz which is 10° higher than magnetic noise 
data [Ref. 16]. The 20 Hz filter cannot completely filter 
out signals of this magnitude. Thus the 4 Hz signal is a 


result of aliasing of the strong ambient 60 Hz signal. 


Pa SLOKES PARAMETERS 
All Stokes parameters demonstrate reasonable agreement. 
That is, Stokes 0 is obviously a? + a* and Stokes 3 is 


Seyiomedy a> - a* as shown in Figures,D.3, Da6, D.13 and D.16. 


47 





I, @eean Site 
The ocean site Stokes 2 and 3 parameters exhibit a 
time dependence. That 1s, in the early a.m. hours the large 
magnitude of Stokes 2 1s the noticeable feature. This 
magnitude decreases at sunrise and through the day until it 
reasserts itself after sunset. A review of Stokes parameters 
2 and 3 over the total observation time (see Technical Report 
NPS-61-83-005) indicates that this falloff is independent of 
frequency. The increase in magnitude of Stokes 2 after 
Sunset is also independent of frequency with the exception 
of some very discrete spectral components. 
2. Land Site 
The land site Stokes 2 and 3 parameters show a time 
dependence with a shorter cycle. That is, the magnitude of 
Stokes 2 1S large near sunset, followed by its decrease 
during the morning hours. The early afternoon hours reflect 
an increase in Stokes 2 until after sunset. Figures 5.1 
through 5.4 show the decreased Magnitude of Stokes parameters 


Beand 5 £or the morning hours. 


C. X AND Y COIL PHASE DIFFERENCE 
1. Ocean Site 
The use of Stokes parameters in the calculation of 
the phase difference dictates that changes in the time 
Series averaged Stokes parameters should be evident in the 


x and y component phase difference. This is the case as 
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shown in Figures D.7 and D.17. The phase difference from 
sunset to sunrise shows a definite decrease with an "eye- 
ball" average of .2 to .4 radians. Daylight hours show 
a difference ranging from .5 to .3 radians with the greatest 
difference near local noon. 
fe sand Site 

The land site phase difference reflects its Stokes 
parameters with a difference of about .3 radians in the 
higher frequencies at all hours and a difference of .1 


radians during the daylight hours. 


ree xX AND Y COIL COHERENCE 
I. Ocean Site 
Time dependence is again a strong factor with high 
degrees of coherence in the ocean Site x and y coils being 
greatest during hours of darkness and diminishing at sunrise. 
Figures 5.5 and 5.6 are examples of this phenomenon. The 
exception to this obServation is at local noon when the 
coherence of the ocean coils is again Very high but diminishes 
Merore andwartcer the noon time data. 
Z. ‘Lande te 
The land site x and y coil coherence reflects the 
previous parameter cycle time dependence. There is a high 
degree of coherence in the lower frequencies at sunrise 
that is followed by diminishing coherence in the higher 


frequencies. Local noon shows an improved coherence at 


oes, 
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Coherence of x and y Coils, Mtry Bay, 18 Aug. SZ, OZ l-O025teL0cal 


Figure 5.5 
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these frequencies but is followed by a degraded coherence 


after sunset. 


Eee CIRCULAR POLARIZATION POWER SPECTRAL DENSITY 
There is no marked contrast between right and left 
@meeular polarization power spectral densities. Both 


calculations provide similar results. 


F, COHERENCE OF CIRCULAR POLARIZATION 

No preference for right versus left circular polarization 
is discernable. The most obvious trend is that of a high 
degree of coherence in the frequency range .05 Hz to .4 Hz 
as shown in Figures 5.7 and 5.8. Generally, coherence of 


the two sites decreases thereafter with frequency. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

In view of the limited data base and the preliminary 
state of the analysis, only tentative conclusions are 
possible at this point. 

ee in the frequency range .05 to .4 Hz, the ocean and 
land data exhibit a high degree of coherence. At higher 
frequencies the coherence rapidly diminishes. 

(2) Both polarized and unpolarized signals were observed 
With a preponderance of high degrees of polarization occurring 
at the lower frequencies. The nature of the polarization 


appears to be a function of both time and frequency. 


B. RECOMMENDATIONS 

The first and obvious recommendation is to collect more 
data. This is already in the planning stage. The second and 
also obvious recommendation is to conduct a more thorough ana- 
lysis of the data on hand. This also is already being planned. 
More concrete recommendations concern data processing. 

One 90-minute real time data tape requires 12 computer 
cylinders of storage space. Thus, for the approximately 
eight hours of data computed and analyzed here, 72 cylinders 
of storage were used. This is too much and is mostly a 
matter of the computer code employed. A little work can 


provide a program that would require less computer storage 


a) 


nf 





space. This will be a real necessity as a data baSe of 
weeks or months is compiled for analysis. It is recommended 
that the present code by "streamlined" to use the IBM 

3033 as a calculation "scratch pad."' Then a few lines of 
"Computer Job Control Language" can be employed to store 

the calculations on a high-density tape for future uSe at 


the Naval Postgraduate School or anywhere in the world. 
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meee JOB (1592,0129),"'FISHER SMC. 1399',CLASS=C 


//*MAIN ORG=NPGVM1. 159 2P, LINES=(60) 


J/*FORMAT PR, DDNAME=PLOT.SYSVECIR, DEST=LOCAL 


feemeon EC FRTXCLGP,PARM. LKED="LISI,MAP,XREF', REGION. GO=2700K 


y7OoXSUTI DD UNIT=SY SDA, SPACE=(CYL, (8,8) ) 


Pei SUT2 Dive oY Sey, SPACE= (CYL, (3,8) } 


J/SYSLIN DD SPACE=(6080,(80,30)),UNIT=SYSDA 


feeee reo SYSIN DD DSN=MSS.SYS3.NONIMNSL.SOURCE(FOURT) , DISD=SHR 


Mf DD.” 

c A LARGE NUYBER OF ARRAYS ARE UPILIZED 

i Withee ine INTENT EON GF MAINTAINING EACH 

oO BEae OF JREIGEINAL DATA OR CALCULATION, 

a peewee TOR RECALL AT ANY TIME FOR FURTHER 

c COMPUTAITIN. EQUIVALENT ARRAYS ARE EMPLOYED 
& PEGBwUSE THE PLOTTING ROUTINE ‘DRAWP' IS 
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NOT ABLE TO HANDLE COMPLEX ARRAYS. 
ARRAYS ‘IPBt,'RTB',"ALABY AND "TITLE! ARE 
USED TO GENERATE THE VERSATEC PLOTTER OUTPUT. 
OCEAN ARRAYS ARE INDICATED VIA THE SUFFIX 
'O', AND LAND ARRAYS ARE INDICATED WITH 'L*. 

COMPIEX*8 XX0(8192) ,YY¥0(3192), 

@s0 (8192) ,co0(8 192) ,'7x9 (3192), 

CRO (8192) ,CTO(8 192) ,CUD (8192), 

CVO (8192) ,CWO(3 192) ,CP0 (3192), 

ZYO (8192) ,COMIL (8192) , COL (8192), 

CUL (8192) ,CVL(8 192) ,CWL(8192), 

COPOL(8192) , XXL (8192) , YYL (8192) , ZXL (8192), 

ZYL (8192) ,CSL(8 192) , BPD (3192), 

BMO (8192) , BPL(8 192) ,BUL(8192), 

SPOL (8192) ,SMOL (8192) ,SP9 (8192), 


Serco 2) ewe (8 192) ,CTL (3192), 
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ei or I 2) 7 eRL (3S 192) , SME (3 19 2) 


Deere NSTON TIME (8192) , FRE) (8192), 


WORK (16384) , FRQ2(8192) 


INTEGER*4Y ITB(12) /12*0/ 


REAL *4 Roe Z2O/2o+0. 0/7, 


CROO (16384) ,CP90(16384) , 2000 (16384), 


COOO (16384) , CT30( 16384) ,2V00 (16384), 


CWOO (16384) ,CTLL(16384), 


CPLL (16384), CRLL(16384) , TOLL (16384), 


CULL (16384) ,CVLL(16384), 


CWLL (16384) ,COMOLL( 16384) , COPOLL (16384), 


SPOO (16384) ,SPLL(16384), 


SMOO (16384) ,SMLL(16 384) 


Beeb ALDAB(22)/° PDZO',*PDYO','*ST10', 


woe O 7 S250,’ STOO', "PHAD', 


ewe ri, PDYL®, "STIL", *ST2L‘, 





'ST3L','STOL',*' PHAL',*SAML!, 


Beemer ere, SUOOT,*SMEL* ,*CPOL', *CMOL'/ 


Pemeers TITLE (12) 


EQUI VALENCE(TITLE(1) ,RIB(5)), (CROO(1), 


mot 1)) , (CP00(1) ,CPO(1)), 


Hemoo( 1) ,coo (1) ) , (CTOO(1) ,CTO(1)) , (CUO0 (1) ,CU0(1) ), 


G70 O( 1) ,CNO (1) ), (CVOO(1) ,CVO(1)), (CILL(1) ,CTL(1}), 


Teer TN Wie), (GR LL( 1 ,pCRL(1)),.(cOLE( 1) ,COLt1} ), 


Pema GU rd) . (CVLL (WW GVL( 1)), (CWLL (1) ,CRL(1)), 


(eer Oli 1).,CO29 (1) ) , (COMOLL (1), 


SOO (iW) ie (SPI9 (1), SP0(1)), 


ety, stl ¢!) ), (>400(1) ,SNO(1)), (SULL(1) ,SML(1) ) 


Dime ao ,YYO,C50,CO0, 2X0, ZYO, CRO, CTO, CUO, CVO, CWO, 


Welt leon leot Lb ,CSi, BPO, BMO,BPY, 


pun, seOL  SMOL,SPO,SPL,SMO, 


CT ln) GRRL, COL ,CVL, COPOL, CNL, 
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CeeCr Ome eMOL/2785928* (0.0,0.0)/ 


DAL AME, FREO, FROZ/28576* 0.0/7 


IFRA ME=8 192 


NR=3 0 


FNR=FLOAT (NR) 


P=0. 0 


DO 70 L1=1,NR 


ie OOmGOIP ENDING WITH STATEMENT 70 ENABLES 


THE PROSRAM TO PROCESS A LARGE AMOUNT OF DATA 


Pen or one ENG THE PROTESS IN BLOCKS.THE DATA 


FOINTS FROM EACH RUN THROUGH THE DO LOOP ARE 


ADDED TOGETHER AND EVENTUALLY AVERAGED BY THE 


NUMBER OF RUNS THROUSH THE DO LOOP. ‘NR! 


Reewoomito DHE NUMBER OF DATA SEQUENCES TO 


BE AVERAGED. 1 SEQUENCE CURRENTLY EQUALS 


8192 DATA POINTS FOR EACH CHANNEL OR 256 


65 





60 


SecONys OF DATA. FOR THES TEST SINE AND 


GGaeuh WAVES OF 4HZ WIL BE USED. 


DO 60 JJ=1,IFRAME 


XXO (JJ) =SIN (P*. 785398163} 


Y YO (JJ) =COS ( P*. 785398163) 


XXL (JJ) =SIN (P*. 785398153) 


YYL (JJ) =COS (P*. 785398163} 


D=P+ 1 


CONTINUE 


ie SP OREONT NG SECTION GENERATES THE TIME 


AND FREQUENCY ARRAYS AND NORMALIZES 


io eneot fcr DATA TO VOLIAGE FORM IN 


PReeeaAl tO WOR THE FAST FOURZTER TRANSFORM 


10 THE FREQUENCY DOMAIN. 


N=8 192 


FN=FLOAT (N) 
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20 


DELRAT=1./32. 


T=FN*DELTAT 


DELT AF=1./T 


DO 20 J=1,N 


TIME (J) =DELTAT* FLOAT (J) 


FREQ (J) =DELTAF* fLOAT (J) 


FRQ2 (J) =ALOG10 ( FREQ (J)) 


CONTINUE 


THE NBMT FOUR STATEMENTS PERFORM 


ANSON THESINPULD TINE SERIES 


DATA. 


CAG Le eOUin (exe , Neal, — 1410..8O0RK) 


CALL FOURT(YYO,N,1, -1,0,W8ORK) 


CALL OUR (XXES NN, 1,-1,0,W ORK) 


CALL, SOU RR(YYL, Njlasi,0, 4 ORK) 


DO 40 K4=1,N 
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40 


XXO (K4) =XXO(K4) /FN 


YYO (K4) =YYO(K4) /FN 


XXL (K4) =XXL(K4) /FN 


YYL (KU) =YYL(K4) /FN 


CONT INUE 


Tee LOOP CALCULATES THE 


OMNGERMALIZED SPECTRAL DENSITIES 


POR = SaNets SITE ORTHISONAL 


COMPONENTS OF THE GEJMAGNETIC FIELD, 


CiageeeiGievipD UAL SITE 2ZROS SPECTRA 


BED WappN COMPONENTS, THE INDIVIDUAL 


CRTHOGONAL COMPONENTS OF RIGHT AND 


LEFT CIRCULARLY POLARIZED SPECTRA, 


PPM WO SETE CROSS SPECTRA YOR RIGHT 


AND GEFT CIRCULAR POLARIZATION, 


pO” 3¢ ITIl=1,N 
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Gee@( UE) =ZXO (IE) + (XXO (II) #CONJG (XXO (II) )) 
ZYO (II) =ZYO (IL) + (YYO (II) *CONJG (YYO(ITI))) 
CSO (II) =CSO(IT) + (XXO (II) *CONJG (YYO (II) )) 
ZXL (IL) =ZXL (IL) + (XXL(II) *CONJG (XXL (IT) )) 
ZYL (II) =ZYL (IL) + (YYL(II) *CONJG (YYL(IT))) 
Gomi =GSL{IL) + (XXL (II) *CONJG (YYL (IT) )) 
BPO (IT) =XXO(II) +((0.,1.) *YYO (IT) ) 

BMO (IL) =XXO(II) -((0.,1.) *YYO (IT)) 
eles ies, sey n(rt)) 

BME (II) =XXL(I1) -((0.,1.} *YYL(ITI)) 

SPOL (II) =SPOL(II) + (BPO (IL) *CONJG (BPL(ITZ))) 
SMOL (IZ) =SMOL(II) +(BMO(II) *CONJG (BML (II))) 
SPO (IL) =SPO(IT) + (BPO (II) *CONJG (BPO (IT) )) 
SPL (IL) =SPh (II) + (BPL(II) *CONJG (BPL(IT) )) 
SMO (II) =SMO (II) +(BMO(II) *CONJG (BMO (IT) )) 


Sb) =Sib(it) +(BML{II) *CONJG (BML (TTI) )) 
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30 CONTINUE 
70 CONTINUE 
THE NEXT LIOP NORMALIZES THE ABOVE 
SPECTRA AND CALCULATES POWER SPECTRAL 
LENS ITIES. 
DO 33 I3=1,N 
ZXO (13) =ZXO(I3) *T/PNR 
ZYO (13) =ZYO(1I3) *T/FNR 
CSO (13) =CSO (I3) *T/FNR 
ZXL (13) =ZXL (13) #0 /FNR 
ZYL (23) =ZYL(I3) *T/FNR 
CSL (13) =CSL(I3) *T/FNR 
SPOL (13) =SPOL(I 3) *T /FNR 
SPO (13) =SPO(I3) *T/FNR 
SPL (13) =SPL(1I3) *T/FNR 


SMO (13) =SMO(I3) *T/FNR 
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Se ( 13) =SML (IS) FT /ENR 


Suet (13) =SHOL(L 3) *T/FNR 


ee CONTINUE 


Pue NEXD LOOP CALCULATES STOKES 


FARAMETERS 


INDIVIDUAL 


COMPONENTS 


THE PLANAR 


FARAMETERS 


DO 44 I4=1,N 


CTitpevcd o OF THe 


SITE ORTHIGONAL 


AND Eie: COUR RENCE OF 


CIRCULAR POLAREZATION 


BolWioen of Les. 


CTO (I4) =(ZXO(I4)+ZYO(I4)) *2./7 


meee) — (2X0 (14 )—-Z2YO (£4) ) ¥2./7(T*CTO(TY) } 


CRO (T4) = (4*CS9 (I4)) / (CII (I4) ¥*T) 


COO (T4) =CSO(I4) /(CSQRT (2X0 (L4) ) *CSORT (ZYO (T4) )) 


CWO (I4) =ATAN2 (A IMAG (COD (I4)),REAL (TIO (T4) )) 


COO (14) =CSQRT (TOO (T4) *CINIG (COO (T4) ) ) 





GTO ( TH) =4 34294 48*CLOG (CTO (14) ) 


CUO (T4) =4.34294 Y8*CLOG(ZX9 (T4) ) 


CVO (I4) =4.34294 48*CLOG(ZYO(I4)) 


CTL (I4) = (ZXL (14) +ZYL (I4)) *2./7 


CPL (I4) = (ZXL (14) -ZYL (14) ) *2./ (T¥CTL (14) ) 


CRL (I4) =(4*CSL(T4)) /(CIL(I4) *T) 


COL (I4) =CSL(I4) /(CSQRT (ZXL (T4) ) *CSQRT (ZYL(T4) ) ) 


CTL (I4) =4. 34294 4U8*CLOG(CTL (I4)) 


CUL (14) =4.3429448*CLOG(ZXL (T4) ) 


CVL (T4) =4.34294 4U8*C LOG (ZYL (T4) ) 


CWL (I4) =ATAN2(AIMAG (COL (I4)),REAL (COL (T4) )) 


COL (I4) =CSORT(COL (L4) *CONJG (COL (T4) ) ) 


COPOL(I4)=SPOL(T4)/ (CSQRT (SPO (T4) ) ¥*CSORT (SPL (TS) ) ) 


COPOL(I4) =CSQORT (COPOL(I4) *CONJG (COPOL (T4) )) 


COMO L(T4) =SMOL(T4)/ (CSQRT (SMO (T4) ) *CSQRT (SML (TY) ) ) 


COMO L(I4)=CSORT (COMOL (T4) *CONJG (COMDL (T4))) 
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4 uy 


SPO (14) =4.3429448*CLOG(SPO(T4) ) 


SPL (I4) =4. 34294 48*C LOG (SPL (I4)) 


SMO (I4) =4.34294 48*CLOG(S49 (T4) ) 


SML (I4) =4.34293 48*CLOG(SML (T4) ) 


COnT INUE 


Vereonroce PLOT OF CALCULATED QUANTITIES 


NPTS=10./DELTAF +1. 


Digsso (SDEETERMINEGS NUMBER.OF POINTS 


HWiemosnnt ©O PLOT THE 0 TO 10 HERTZ RANGE. 


ITB (2) =0 


ITB ( 3) =20 


ITB (4) =10 


ITB (6) =1 


ITB (7) =1 


ITB ( 12) =1 


RTB ( 1) =0 





RTB ( 2) =0 


RTB (3) =ALAB(1) 


READ (5, 3000) TITLE 


CALL DRAWP(NPIS,FRQ2,CU0),ITB,RTB) 


ITB (6) =1 


RTB (3) =ALAB (2) 


READ (5,3000) TITLE 


CALL DRAWD(NPTS ,FRO2,CVOD,ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB (3) 


READ (5,3000) TITLE 


Sper nnwe (els , © REOQ,C POI, TTS, RIB) 


ITB (6) =1 


RTB (3) =ALAB (4) 


READ (5, 3000) TITLE 


CALL DRAWP(NP?DS ,FREQ,CROD,ITB, RTB) 
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ITB (6) =1 


RTB (3) =ALAB(S5) 


READ (5, 3000) TITLE 


CALL DRAWP(NPTS,FREQ,CROO(2) , ITB, RIB) 


ITB (6) =1 


RTB (3) =ALAB(6) 


READ (95,3000) TITLE 


ih URAWP(NPTS , FROZ,CIO), [TB, RTB} 


ITB (6) =1 


RTB (3) =ALAB(7) 


READ (5, 3000) TITLE 


CALL DRAWP(NPIS ,FREQ,CWOO0,ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB (8) 


Ren (5,3000) TITLE 


CALL DRAWP(NPIS, FREQ,COO),ITB,RTB) 


figs: 





ITB (6) =1 


RTB (3) =ALAB(9) 


READ (5,3000) TITLE 


CALL DRAWP(NPIS, FRQ2,CULL,ITB,RTB) 


ITB (6) =1 


RTB (3) =A LAB (10) 


READ (5, 3000) TITLE 


Gaul DRAWP(NPIS ,FROZ,CVLL,ITS, RTS) 


ITB (6) =1 


RTS (3) =ALAB(11) 


READ (5,3000) TITLE 


Grol DRARP(NETS , FREQ, CPLL,I I TSB,RTB) 


ITB (6) =1 


RTB ( 3) =ALAB(12) 


READ (5, 3000) TITLE 


Chie DRMWPENPTS ,FREQ,CRLL,1TB,RTB) 
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ITB (6) =1 


RTB (3) =ALAB(13) 


READ (5, 3000) TITLE 


Sab PRAWP{NPTS ,FREQ,CRLL (2) ,iITB, RIB) 


ITB (6) =1 


RTB (3) =ALAB(14) 


Pare, 3000) TITLE 


Sale DRAWP(NPIS, FRQO2Z,CILL, 1[TB, RTB) 


ITB (6) =1 


RTB ( 3) =ALAB (15) 


READ (5, 3000) TITLE 


CALL DRAWP(NPTIS ,FREQ,CWLL,ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB(16) 


READ (5, 3000) TITLE 


CALL DRAWP(NPIS, FREQ,COLL,ITB, RTB) 
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ITB (6) =1 


RTB (3) =ALAB(17} 


READ (5,3000) TITLE 


Sri DRAWP(NPIS,FRO2,SPO),1TB,RTB) 


ITB (6) =1 


RTB (3) =ALAB(18) 


READ (>, 3000) TITLE 


Orr DRAWP (NETS ,FROZ2,SPLL,ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB( 19) 


READ (5, 3000) TITLE 


Cape DRAWE{NPIS, PROZ,SM0O), i1TB, RTB) 


ITB (6) =1 


RTB (3) =ALAB (29) 


Reno, 2000) TITLE 


Sere OR AWE (NPIS ,FROZ,SMLL,ITSB,RTB) 
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ITB (6) =1 


RTB (3) =ALAB(21) 


READ (5,3000) TITLE 


Brel DRAWP(NOTS ,FREQ,COPOLL,IT B&B, RTB) 


ITB (6) =1 


RTB ( 3) =ALAB(22) 


READ (5, 3000) TITLE 


mr DRAWP(NPIS, FREQ, CONILL,ITB5,RTB) 


3000 FORMAT (6A8) 


5.0 P 


END 


/* 


ByeOs.5YSIN DD * 


meet POWER SPECTRAL DENSITY OF SIN 4AZ 


IN XX0,30 AVGS, (REF DB VS LOG FREQ) 


oo POWER SPECTRAL DENSITY OF [OS 4HZ 
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fee, 30 RVGS, (REF DB VS LOG FREQ) 


mot STOKES I/STOKES J OF SIN WAZ 


IN XXO AND COS 4HZ IN 


YYO, 30 AWSS 


Meee STORES 2/STOKES 0 OF SIN 4HZ 


IN XxXO AND COS 4HZ IN 


YY¥O;20SAVsS 


meer STORES 3/7STOKES 0 OF SIN 4HZ 


IN XXO AND COS 4HZ IN 


YYO,30 AVSS 


Mier STOKES O OF SIN 4HZ IN XXO AND 


COS 4HZ IN YYO,30 AVSS 


mol PHASE OF SIN 4HZ2 


AND COS 4HZ IN YYO, 30 


Mest COHER OF SIN 4HZ 


AND COS 4HZ IN YYO, 30 


S(REP DB VS LOG) 


IN XxXO 


AVGS 


IN XXO 


AVGS 


Tie POWER eSPECTRAL DENSITY OF SIN 4H2Z 


IN XXL,30 AVGS, (REF DB VS LOG FREQ) 


Me>t POWER SPECTRAL DENSITY OF [COS 4HZ 
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mor yn,30 AVGS. (REF DB VS LOG FREQ) 


MEST STOKES 1/STOKES 9 OF SIN 4AZ 


IN XXL AND COS 4HZ IN YYL,30 AVSS 


meee STORES 2/STOKES J OF SIN 4HZ 


MP seL AND COS 4HZ IN YYL,30 AVSS 


Moot STOKES 3/STOKES 0 OF SIN 4HZ 


Mee xX AND COS 4SHZ IN YYL,30 AVSS 


Mise STOKES 0 OF SIN 4HZ IN XXL AND 


S@oethHzZ IN YYL,30 AVSS, (REF DB ¥S LOG) 


Mest PHASE OF SIN 4HZ IN XXL 


AND COS 4HZ IN YYL, 30 AVGS 


Meee COMER CF SIN 4HZ IN XXL 


Pee COS 4HZ IN YYL, 30 AVGS 


meet RE CLiRC POLARIZATION PSD OF SIN GHZ IN 


mo BND COS SHZ IN YY9, (DB VS LOG),30 AVSS 


mint CLERC POLARIZATION PSD OF SIN 4HZ IN 


on 





joo AND COS 4HZ IN Y¥L, (DB VS LOG),30 AVGS 


meoe LEFT CIRC POLARIZATION PSD OF SIN 4HZ 


IN XXO AND COS 4HZ IN YYO (DB VS LOG) ,30 AVGS 


Meee LEFT CIRC POLARIZATION PSD OF SIN 4HZ 


IN XXL AND COS 4HZ IN YYL (DB VS LOG),30 AVGS 


Meod COHER RE CIRC POLARIZATION OF SIN 4HZ 


mien kO,XXL AND COS 4HZ IN YYO,YYL, 30 AVGS 


meer COHER LEFT CIRC POLARIZATION OF SIN 4HZ 


IN XXO,XXL AND COS 4HZ IN YYO,YYL, 30 AVGS 


/* 


// 
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FREQUENCY (HZ) 


Test Ocean Stokes Parameter Three 


Figure A.6 
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FREQUENCY (HZ) 


Test Land Stokes Parameter Three 


Figure A.15 
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FREQUENCY (HZ) 


Test Land Coherence of x and y Coils 


Figure A.17 
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NASS STORAGE PRISKAU 


feeeoime og CB (1592,0129)," FISHER SMC.1399',CLASS=F 


//*MAIN ORG=NPGVM1. 15922 


// EXEC FORTXCLG 


Peer te sitslN DD DSN=4SS.SYS3.NONIMSL. SOURCE (FOURT) , DISP=SHR 


// DD * 

c THE ARRAY 'IN' WILL BE USED TO 
c RECEIVE THE DATA PASSED FROM 

C THE SUBROUTINE 'READ' AND THEN 
C TRANSFERRED TO THE APPROPRIATE 
C XXX OR YYY ARRAY. 


INTEGER*2 IN(15) 


COMPLEX*8 XXX (8192), YYY (8192) 


DIMENSION FREQ(8 192) ,TIME(8192) , WORK (16384) 


DATA XXX,YYY/153 84*(0.0,0.9) / 


DATA TIME, FREQ/16384*Q.0/ 
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Se, 


Ms rPOLLOWENG SECTION READS THE 


Perot onk SPGONDS JF COMPUTER TAPE 


AND DISCARDS THIS DATfA. 


Pmi= 1°92 


DO 55 JJ=1,I1TL 


eee RD (20,2N,200, TREC, IRR) 


CONTI NUE 


IFRAM F=8192 


NR=20 


FNR=F LOAT (NR) 


DO 70 L1=1,NR 


ome NEAT LOOPeRERDS THE COMPUTER 


Par etsrNs THEMPROSRAM PROVIDED 


Breen. Tit —SRaeON OF NAVAL POST- 


GRADUATE SCHOOL. 


DO 60 JJ=1,IFRAME 


LO 





Chere ( 20,6 71000 ,T REC, IRR) 


Kee gs)l= = Ni(2} 


YYY(J J) =IN (3) 


CONTI NUE 


N=8192 


FN=PLOAT (N) 


DELTA T=". /32. 


T=FN* DELTAT 


Dewi B= 12 /T 


DO 20 J=1,N 


TIME (J) =DELTAT*F LOA (J) 


FREQ (J) =DELTAF*F LOAT (J) 


THe oNExE 4 STEPS CONVERT THe 


Deite tO VILAGE ANDSENSURES THAT 


NO ERRONEOUS DATA HAS BEEN INTRODUCED 


Oe THEA RAYS < 
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XXX (J )= (XXX (J) -2 045.5) *10. 2045.5 


XXX (J)=REAL (XXX (J) ) 


YYY (J)= (YYY(J)-2 045.5) *10. 2045.5 


YYY(J)=REAL (YYY(J)) 


20 CONTINUE 


CALL FOURT(XXX,N ,1,71,0,#ORK) 


Gwe cOURT{(YYY,N ,1,71,0, WRK) 


DO 40 K4=1,N 


XXX (K 4) =XXX(K4) / FN 


YYY(K4) =YYY(K4) /FN 


GO CONTINUE 


ida NEKT LOOP APBLIES THE SYSTSa 


feeitoein S£UNCTRON TD THE TRANSFORMED 


FREQUENCY DOMAIN DATA. THE TRANSFER 


PONCTTON CONVERTSMVOPTS TO NANOTESLAS. 


DO 9 L=1,N 


im"? 





FRQ= F FEQ (L) 


TP (ae. LE. 152) so TO 2 


KXX(L)=XXX (L) /(105.5-3. 14*PRQ) 


meaeEy—Y iy (lL) 7181.3 2-7.588*F RQ) 


GO TO 8 


meme. LEL10.)S) To 3 


XXX (L)=XXX (L) /(5 .958*FROQ-30. 97) 


meee) —Y YY (L) /(7 .165*FRI-39.99) 


GOP sO S 


[P(e ets oy) oo Lo 4 


exe L =k WLP(3 292 = FRO-5 . 31) 


eer = erere(L) 78 2252 O- 10.85) 


LE MGR Coes GO TO"5 


KCL ) =X XR WL) / (2 .6311*FRO+0. 14667) 


emer =m (yey (3 200M PRO-1.55) 





ROe Less) CO TOG 


XXX (L)=XXX (L) /(2 .6311*FRO+0. 14667) 


meer (L)=Y YY (L) /(2.702*FRQ) 


GO TO 8 


MXX(L)=XKX(L) / (2.7 2* FRO) 


GO TO 7 


CONTINUE 


CON NU E 


iheoeeet WRITE STADEMAENTS SEND 


tne  CONVERL ED DATA TI MSS 


FOR FUTURE MANIPULATION AND RECALL. 


Reo ct oS hl OF STRTERENTS ALSO 


PROVIDES THE USER A DIAGNOSTIC. 


WRITE (21) XXX 


WRITE (6,*) XXX (1) , XXX (8192) 
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WRITE (21) YYY 


Woieeemtoyrr ye erecip , YYY(S 719 2) 


70 CONTINUE 


ENDFILE 27 


SUR OY 


END 


Seenecul iNe BO(LUN, 10, 1RS, TREC, IRQ} 


fio nee DURE FURNISHED BY MR. FIM STANTON, 


DEPARTMENT OF OCEANOGRAFHY. 


READ DATA FROM IUN, ALLIGN , CHECK & RETO RN 


IUN=TAPE NUMBER, £G 2) 


IO=INTEGER*2 ARRAY, 16 LONG, 


(VALUES 0-4095, SUBTRACI 2048) *5/2028. GIVES VOLTAGE 


110 





PRS= NUMBER OF RESINCS ALLOWED (ERRORS) 


ITREC= COUNTER OF RECORDS (FRAMES OF DATA) 


BLeGh 2 BITS, Sess == XECORD 


800 BPI TAPE UNLABLED 


meo- NUMBER OF ACTUAL RESINCS (ERRORS) 


INTEGER * 2 I9(16),1IP(16) 


DATA IRR /0/ 


Pe tREeSEO.0) IS=0 


FORMAT (16A2) 


me (Ss NEsO) Sd TO 50 


READ (IUN,20,END=900) IP 


TREC=IREC+1 


io eot | 


2. 





50 


ays) 


70 


fects. Le. i?) GO TO 50 


READ (I0N,20,END=900) IP 


TS=1 


TRES= TREC+1 


a= 0= MASK (IP (IS) ,3,9)+1 


Meee (6,55) L[CH,1IS,IUN, TREC 


Perm (" KRESYNCING ICH,IS,IUN,IREC ',418) 


im (SCHONES1) SO TO 40 


DO 100 I=1,16 


GGWL)=LSHiPece (15) -4) 


m@m=i ASK (2P (1S) ,3,0) +1 


foaeredsEOQ<ijpeso TO 80 


TE R=IER+1 


Wiese (6,70) LUN,IREC ,I,1CH,IER 


FORMAT (*UNIT', I3,"RECORD',1I6,'CHAN & DATA CH 


2 


',214, 





&0 


100 


"ERRORS ',I7) 

IS=1IS+1 

ers bDeiapes0. TO. 109 
READ (IUN,20,END=900) IP 
IS=1 

IREC=IREC+1 


CONTINUE 


De TERS EOLOy so Tr 7150 

IRR=IRR+1 

Me IRR-LT.1RS) GO TO 129 

WRITE (6,110) 

FORMAT ('1 STOPPED IN SUB RD BECAUSE 
O-MLGR AG Tats, OAT. L170!) 

TRQ=IRR 


Soe 





120 


130 


150 


900 


S110 


CONT INUE 


WRITE (6, 


130) IREC, IRR 


eA eG unbovvomnt FRAME’ ,T6,*WITH TOTAL ERRORS', [7) 


TRQ=IRR 


eee TO 50 


CONTINUE 


RETU &N 


WRITE (6, 


910) IUN,I REC 


FORMAT ('1 END OF UNIT ',1I3,' AT REC ',17) 


> Ose 


END 


FUNC TION 


ISHIFT (IN, NPLC) 


RETURNS SHIFTED VALUE OF I*2 WORD IN 


Ji foe RIGHT SHIFT 





INTEGER * 2 IN 


ime— = N 


foetrs bre) LP=i Pt65536 


IF (NPLC.LT.0) GO TO 30 


ISHIFT=IP/(2**I ABS(NPLC) } 


RET U RN 


ISHIFT=IP*(2**I ABS(NPLC) ) 


TF (ISHIFT.GT.65535) ISHIFT=MOD (ISHIFT, 65536) 


REI U RN 


END 


FUNCTION INASK (IN, IBL, IBR) 


Leia lLeeanwORD IN OUTSIDE BITS LBL & IBR 


INTEGER * 2 IN, IO 


TO=IN 


Git tote O20) 930 T) 50 


SEIU 





t= Psnirt (iN,18 R) 


LO=I1T 


50 Pe-seshir T(iO,LBL~15-IBR) 


nO= 1 .P 


IMAS K=ISHIFT (I) ,15-IBL) 


RETURN 


END 


f* 


V7oceer 2 IFO001 DD UNIT=3330V,MSVSP=PUB4A, DISP=(NEW,CATLS) , 


v7, DSN=Wee. SIS92.GMDIIA, 
7 / DCB= (REC FM=VBS, BLXSIZE=4096, LRECL=4092) , 
/ SERCE=(c YL,(8,8) ) 


//GO.FT2Z20F001 DD UNIT=3400-4,VOL=SER=GMDT1A,DISP=(OLD,PASS), 


V/s ASEL—=(1,NL, ,LN) , 
u/ DCB= (REC FM=FB,LRECL=32, BLKSIZ£=512, DEN=2) 
af 


LG 





feeeen iis JOR (1592,0129),'FISHER SMC. 1399',CLASS=F 


//*H8AIN ORG=NPGVM1. 159 2P,LINES=(60) 


//*FORMAT PR,DDNAME=PLOT.SYSVECTR, DEST=LOCAL 


Pyeeerese FREXCLGP,PARM. LKED="LIST,MAP, XREF',R EGION.GO=270 0K 


7s SUT DD UNIT=SY SDA, SPACE=(CYL, (8, 8) ) 


yo SUT2 DUE —o i 5D A, SPACE=(CYL, (8,8) ) 


J/STUSLIN DD SPACE=(6080, (80,39)),UNIT=SYSDA 


wae oto N DD DSN=4SS.SYS3.NONIMSL.~SOURCE (FOURT) , DISP=SHR 


// 1 

C Pitney oO sScnleirons OF THE 

cS VARIOUS STEPS AND LOOPS ARE 

c CONTeaieo EN APPENDICES A AND 8B. 


INTEGER*2 IN (15) 


GUumemrsatS KKOCS 192) ,YYO(8 192), 


C CSO (8192) ,CO0(8 192) ,2X9(8192), 


de yy 





CRO (8192) ,CTO(8 192) ,CU0 (3192), 


CVO (8192) ,CWO(8 192) , CPD (8192), 


ZYO (8192) ,COMOL (8192) ,COL (8192), 


SMS 12) PCV Li S192) CNL (3192), 


COPOL(8192) , XXX (8192) , YYY (8192), 


XXL (8192) ,YYL(8192) ,2XL (3192), 


ZYL (8192) ,CSL(3 192) ,BP9 (8192), 


BMO (8192) ,BPL(8 192) ,BML(3192), 


Erper (atS2),SMOL (8192) ,S22 (8192), 


ere e 19 2)e7SN0(38 192) , CTE (8192), 


@Pm¢ 8192) -CRL (9192) , SMBS 19 2) 


DIMENSION TIME(8192) , FRED (8192), 


WORK (16384) , FR) 2(8192) 


INTEGER*4 ITB(12) /12*0/ 


REAL *4 moet 2 8j728*0.07,TR00 (16384), 


CPOO (16384) ,CU90(16384) , 


ies 





COOO (16384) ,CTI0(16 384) ,2VO0 (16384), 


CWOO (16384) , CILL(16384) , 


CPLL (16384) ,CRLL(15 384) ,20LL(16384), 


CULL (16384), CVLL(16384), 


CWLL (16384) , COMOLL( 16384) , COPOLL (16384), 


SPOO (16384) ,SPLL(15384), 


SMOO (16384), SMLL(16384) 


Beet, ALAB (22) 7* PDXO! ,*PDYO’,*ST10°, 


woe" | * ST30',* STOO! ,* PHAD', 


Maro", "PDxL®,"PDYL", "STIL',*ST2L', 


women | SOL", * PHAL',*SAML*, 


Po eee on, SUOO®, “SILL',*CPOL',*CMOL* / 


REAL*8 TITLE (12) 


Retr yAbmNCE (TITLE (1) ,RITB(5)), (CROO(1), 


Geet) ) , (e200 (1) ,CeO(1)), 


(G0G 0111) Geo (1) ), (<1Te0(1) ,CTO(1)), (CU00(1),CU0(1)), 


i ie, 


| 





Q 


fenowr ty ,enu (1) ), (£V00(1) ,5V0(1)),(CTLL(1) ,CTL(1)), 


Ceo, eri (i) ), (CRLL(1) ,CRL{(1)), (COLL (71) ,COL(1)), 


(CULL(1) ,CUL (1) ), (CVLL(1) ,CVL(1)), (CWLL(1) ,CWL(1)), 


(eer cLi{1),CGP2L(1)), (COMOLL(1), 


Somgi(l)),(SPI9 (1),SP0(1)), 


eerie or et), (oOo (i), sO (1)),(SHLL(1) ,SML(1} ) 


DATA XX0O,YYO,CS0,CO0,2X0, ZY0O, CRO, 


wre cu, CVO, CH) , XXX, YYY, 


eaevil 7s, ZYL,CSL,BPO,BMO,BPL, 


BML, SPOL,SMOL, SPO,SPL,SM), 


wel u,enk, CoOL ,CVL,COPOL,CWL,COPOL, 


Come L7294912*(0.0,9.0) / 


Dia lier noo, FROZ/S2ZU575*0. 0/ 


£PL= 192 


DO 55 JJ=1,ITL 


CALL RD(20,IN,200,IREC,IBR) 


E20 





CONTINUE 

I FRA ME=8 192 

NR=20 

FNR= FLOAT (NR) 
THE FIRST STATEMENTS OF LOOP 70 
RECALLS DATA PREVIDUSLY STORED IN 
THE IBM 3033 MASS STORAGE SYSTEM. 

DO 70 L1=1,NR 

READ (21) XXX 

READ (21) YYY 

DO 43 II=1,IFRAME 

XXL (IT) =XXX (II) 

YYL (IZ) =YYY (II) 

CONTINUE 

DO 60 JJ=1, IFRAME 


GREL a Da Z07e N00, TREC,IRR) 


sa 





60 


XXO (JJ) =IN(2) 


YYO (JJ) =IN(3) 


CONT INGE 


N=8 192 


FN=F LOAT (N) 


DELTAT=1./32. 


T=FN*DELTAT 


DELTAF=1./T 


DO 20 J=1,N 


TIME (J) =DELTAT* FLOAT (J) 


FREQ (J) =DELTAFP* FLOAT (J) 


FRQ2 (J) =ALOG 10 ( FREQ (J) ) 


XXO (J) = (XXO (J) - 2045.5) *5./2045.5 


XXO (J) =REAL (XXO (J)) 


YYO (J) = (YYO(J) - 2045.5) *5. 72045.5 


YYO (J) =REAL(YYO (J)) 


HP 


i 





20 


40 


CONTINUE 


CALL FOURT(XX0O,N,1,-1,0,WORK) 


CALL FOURT({YYO, N,1,-1,9,WORK) 


DO 40 K4=1,N 


XXO (KU) =XXO(K4) /FN 


YYO (K4) =YYO(K4) /FN 


CONTINUE 


DO 9 L=1,N 


FRO= FREQ (L) 


Mme sO. bES15.)S0 TO 2 


XXO (L) =XXO(L)/(-81. 36 +17.066*FRO - 0.5389*FROQ*FRQ) 


YYO (L) =YYO(L) /(-63. 435 +11.7106*FROQ-0.23279*F RO*F RO) 


GO TO 8 


TRON .10.)50 T) 3 


XXO (L) =XXO(L)/(-18.486 + 4.9029*FRQ) 


Yer 1.) = PO) 7 (6 ser8 35 +71. 7106*PRO -0.23279*FRQ* FRQ) 


5 





GOuT.0 8 


meee RO. LE./.5)S0 TO 4 


XXO (L) =X XO (L) /(3- 94 7*FRQ-9. 368) 


YYO (L) =YYO(L) /( 4. 295*FRQ-11.837) 


Eon. O 8 


mug rOs LEgo.)sd TO 5 


XXO (Ll) =XXO(L)/(2.8105¥*FRQ-0.5657) 


Pave) = Y YOUL) / (2.80 UF FRQ-.50) 


GCOnT.0 8 


XXO (L) =XXO(L) /( 2. 723*FRQ) 


WeaorcL) =YYO{L)/ (2.71 7*FRQ) 


GO TO 8 


CONT INOE 


CONTINUE 


DO 30 II=1,N 


ZXO (IZ) =ZXO (II) + (XXO (II) *CONJG (XXO (II) )) 
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30 


Aoeweiijy=2rO( rE) t(yYYO(LL) *CONJG (YYO (ITI) )) 


CSO (II) =CSO (II) +(XXO (II) *CONJG (YYO (ITI) )) 


Peel P= GRE (ht) + (XX L(LL)“*COINJG (XXL (IT) )) 


ZYL (II) =ZYL(IL) + (YYL(IL) *CONJG (YYL(ITI))) 


CSL (II) =CSL(ILT) + (XX L( II) *CONJG (YYL(II))) 


BPO (II) =XXO(II) +((0.,1.) *YYO (IT) ) 


mec LL) =XXO(II) -((0.,1.) ®YYO (TI) ) 


biieg@ut) =XXiL(it) +((0.,1.) *YYL(i1I) ) 


Bae (IL) =XxXL(IiL)-((0.,1.-) *YYL (ZT) ) 


SPOL (II) =SPOL(I 1) + (BPO(1I) =CONJG (BPL(1ITI))) 


Siete) —SuOL(i I) + (BNO (TL) *CONJG (BML (IT) ) ) 


eee) — SFO (LL) + (BPO (LL) *CONJG (BPO (ITI) )) 


Sete) = SPE(IL) + (BP L(il) FCONJG (BPL(I1I) )) 


SMO (II) =SMO(II) + (BMO (II) *CONGG (BMO(ITI))) 


SML (IT) =SML(IT) + (BUL (IL) *CONJG (BML (ITI) ) ) 


CONTINUE 





70 


23 


CONTINUE 


DO 33 I3=1,N 


ZXO (13) =ZX0(1I3) *T/FNR 


ZYO (13) =ZYO (13) *T/FNR 


CSO (I3) =CSO(I3) *T/FNR 


ZXL (13) =ZXL(I3) *T/FNR 


ZYL (13) =ZYL(1I3) *T/FNR 


CSL (13) =CSL(I3) *T/FNR 


SPOL (I3) =SPOL(I 3) *T/FNR 


SPO (I3) =SPO(I3) *T/FNR 


SPL (13) =SPL(1I3) *T/FNR 


SMO (I3) =SMO (13) *T/FNR 


SML (13) =SML(1I3) *T/FNR 


SMOL (13) =SMOL(I 3) *T/FNR 


CONTINUE 


DO 44 I4=1,N 
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CTO (I4) = (ZXO (14) +ZYO (14) ) *2./T 


CPO (IU) = (ZXO (I4) -ZYO(I4)) *2./ (T*CTO(I4)) 


CRO (14) = (4*CS9 (I4)) / (CTO (I4) *T) 


COO (14) =CSO(I4) /(CSQRT (ZX) (T4) ) *CSQORT (ZYO (T4) )) 


CWO (I4) =ATAN2(AIMAG (COO (IT4)) , REAL (COO (T4))) 


COO ( T4) =CSQRT (COO (T4) *CONIG (COO (T4) ) ) 


CTO (T4) =4.34294 48*O LOG (CTO (T4) ) 


CUO (14) =4.34294 48*C LOG (ZXO (I4) ) 


CVO (14) =4.34294 U8*CLOG(ZY9(T4)) 


CTL (14) = (ZXL (14) +ZYL(I4)) *2./T 


CPL (I4) = (ZXL (14) -ZYL(1T4)) *2./ (T*CTL(T4) ) 


CRL (IU) = (4*CSL(I4)) /(CIL(I4) *T) 


COL (14) =CSL(I4) /(CSQRT(ZXL(I4) ) *CSQRT (ZYL (T4) )) 


CTL (14) =4.3429448*C LOG (CIL (I4) ) 


CUL ( 14) =4.34294 U8*CLOG(ZXL (I4) ) 


CVL (I4) =4.34294 48*CLOG(ZYL (T4) ) 


a2 7 





Gy 


CWL(I4) =ATAN2(A IMAG (COL (T4)), REAL (COL (14) )) 


COL (I4) =CSQRT (COL (14) *CONJG (COL (I4)) ) 


COPOL(I4&) =SPOL(I4)/ (CSQORT (SPO (T4) ) *CSQORT (SPL(T4) ) ) 


COPOL(T4&) =CSQORT (COPOL (14) *CONJG (COPOL (I4))) 


COMOL(I4) =SMOL(I4U)/(CSIRT (SMO (I4)) *2SORT (SML(I4)) ) 


COMOL(I4&) =CSORT (COMOL(I4) *CONJG (COMOL (T&) ) ) 


SPO (T4) =4.34294 48*CLOG(SP9 (T4) ) 


SPL (14) =4.34294 YU8*T LOG (SPL (I&4) ) 


SMO (T4) =4.34294 48*C LOG (SMO (I4)) 


SML (14) =4. 34294 4U8*CLOG(SML(TI&) ) 


CONTINUE 


Nels —1Ue—e/DELTAY +1. 


ITB ( 2) =0 


ITB (3) =20 


ITB (4) =10 


ITB (6) =1 


128 





ITB (7) =1 


ITB (12) =1 


RTB ( 1) =0 


RTB (2) =0 


RTB (3) =ALAB(1) 


Fem (>, 3000) TITLE 


CALL DRAWP(NPTS ,FRQ2,CUO2,ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB(2) 


Meee, 2000) TIT LE 


CAGE DRAWE(N PES ,FROZ2,CYO),1TB, RTB) 


ITB (6)=1 


RTB (3) =ALAB (3) 


READ, 3000) TET LE 


Cac UrmAee (NEES ,YREO,CPO), 173, RTB) 


ITB (6) =1 


L208, 





RTB (3) =ALAB(4) 


READ (5,3000) TITLE 


CALL DRAWP(NPIS,FREQ,CRO),ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB(5) 


READ (5, 3000) TITLE 


LL, DRAWN PCNPrs ,PREO,CRO) (2) , ITB, RTS) 


ITB (6) =1 


RTB (2) =ALAB (6) 


meno (>, 3000) TISLE 


CALL DRAWP(NPTS ,FRO2,CTO),iTB, xTB) 


ITB (6)=1 


RTB (3) =ALAB(7) 


Bend (5,3000) TITLE 


Cae De awe(NPTS > Fre0,chwO), LTB,RTB) 


ITB (6) =1 


1S 0 





RTB (3) =ALAB(8) 


READ (5, 3000) TITLE 


CALL DRAWP(NPIS,FREQ,CO0O0,ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB(9) 


READ (5, 3000) TITLE 


Ger DRAWP(NPTS ,FRO2,CULL, ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB(10} 


READ (5, 3000) TITLE 


Gres DRAWE (NPIS, FROZ,CVLL,ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB(11) 


READ (5, 3000) TITLE 


CALL DRAWP(NPTS ,FREQ,CPLL,ITB, KTB) 


ITB (6) =1 


oe 





RTB (3) =ALAB(12) 


READ (5, 3000) TITLE 


CALL DRAWP(NPTS, FREQ, CRLL, ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB(13) 


Ree (5, 3000) TETLE 


CALL DRAWP(NPIS, FREQ,CRLL(2) , ITB, RTB) 


ITB (6) =1 


RTB ( 3) =ALAB(14) 


ReaD (5, 3000) TITLE 


CALL DRAWP(NPIS,FRQ2,CTLL,ITB, RTB) 


ITB (6) =71 


RTB (3) =ALAB(15) 


ReaD (6, 3000) TITLE 


Gail DRAW (NPS ,PREO, CULL, ITB, RTB) 


ITB (6) =1 


N74 





RTB (3) = ALAB( 16) 

READ (5, 3000) TITLE 

CALL DRAWP(NPTS,FREQ,COLL,ITB, RTB) 
ITB (6) =1 

RTB (3) =ALAB(17) 

READ (5, 3000) TITLE 

CALL DRAWP(NPTS ,FRO2,SP09,ITB, RTB) 
ITB (6) =1 

RUB (3) =ALAB(18) 

READ (5, 3000) TITLE 

CALL DRAWP(NPIS,FRQ2,SPLL,ITB, RTB) 
ITB (6) =1 

RTB (3) =ALAB(19) 

READ (5, 3000) TITLE 

CALL DRAWP (NPIS ,FRQ2,SM0),ITB, RTB) 


ITB (6) =1 





RTB (3) =ALAB(20) 


READ (5, 3000) TITLE 


Sele OURAWE(UPES , FROZ,SMLL,ITS, RTB) 


ITB (6) =1 


RTB (3) =ALAB(21) 


READ (5, 3000) TITLE 


CA DRA weir lS PEREQ, COPILL, ITB, RTB) 


ITB (6) =1 


RTB (3) =ALAB(22) 


READ (5,3000) TITLE 


CALL DRAWP(NPTS , FREO,COMOLL,ITB, RTB) 


3000 FORMAT (6A8) 


STOP 


END 


Sueneud muse RD(LUN, LO, ERS, IREC, IRQ) 
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PoeoeesOCEDURE FURNISHED BY MR. [IM STANTON, 


DcPARTMENT OF OCEANOGRAPHY. 


READ DATA FROM JUN, ALLIGN , CHECK & RETURN 


IUN=TAPE NUMBER, EG 2) 


IO=INTEGER*2 ARRAY, 16 LONG, 


(WREUES O0-4095 , SUBTRACT 2048) *5/2028. GIVES VOLTAGE 


IRS= NOMBER OF RESINCS ALLOWED (ERRORS) 


TREC= COUNTER OF RECORDS (FRAMES OF DATA) 


SEOGK 512. B815TS, 32 BITS = RECORD 


800 BPI TAPE UNLABLED 


Pe ieee or ACTUAL RESINCS (ERRORS) 


INTEGER * 2 I10(16),IP(16) 





20 


40 


50 


55 


DATA IRR /0/ 


IF (TREC.5£Q.0) IS=0 


I ER=0 


FORMAT (16A2) 


IF (IS.NE.0) G2 TO 50 


READ (IUN,20,END=900) IP 


PREC=IREC+1 


ITS=I S+1 


Peeeielos Gt.17) sO TO 50 


READ (IUN,20,=ND=900) IP 


IS=1 


IREC=IREC+1 


TGRe I MASK (IP (15) ,3,0) +1 


WeeereGorm o>) I[H,IS ,LUN,IREC 


Rotate (* RESYNCING ICH, 1IS,i1UN,IREC 


a 9S) 





Pp mCLehea ret). SO: TO-40 
DO 100 I=1,16 
IO(I)=ISHIFT (IP (IS) , 4) 
ICH=I MASK (IP (IS) ,3,0) +1 
MaMe@mGH.—EQ.1) GD TO 80 
IER=I ER+1 
WRITE (6,70) IUN,IREC ,I,ICH,IER 
70 FORMAT (‘UNIT',1I3,'RECORD',1I6,'CHAN & DATA CH ',2I4, 
"ERRORS ',17) 
80 IS=IS+1 
fos. UTP SO TO 109 
READ (IUN,20,END=900) IP 
IS=1 
IREC=IREC+1 


100 CONTINUE 





IF (IER.E0.0) GO TO 150 


IRR=IRR+1 


Mia (RR. LL. iRS) GO TO 120 


WRITE (6,110) 


1106 newer (* beSTOPPED IN SUB RD BECAUSE 


moral RR. Gln to, Al Eto") 


TRQ=IRR 


S10 P 


120 CONTINUE 


WRITE (6,130) IREC, IRR 


130 Bonne heot NGC AT FRAME’ ,16,"*WITH TOTAL ERRORS',I7) 


TER=0 


IRQ=IRR 


Gemiro 5 0 


150 CONTINUE 


RETU EN 


ise 


fi 


u 





a0 


oO 


30 


WRITE (6,910) IUN,IREC 


BORMAT ('T Pippeersunee *,t3,° AT REC *,17) 


SOP 


END 


MUNG TLON Poitier l (IN, NFLC) 


POUURioescH ETE VALUE OF I*2 WORD IN 


iMoney Ristit Sai FT 


Pi CER. s 22. CN 


Ip=iIN 


IF (IP.LT.0) IP=IP+ 65536 


MAPLE C.LT.0) GO TO 30 


heh. P/(2-FEABS (NPLC) | 


Role RN 


ISHIFT=1P*(2**f ABS{NPLC) ) 


IF (ISHIFT.GT.65535) ISHIFT=MOD (ISHIFT, 65536) 





5 0 


/* 


RET U RN 


END 


FUNCTION IMASK (IN, IBL, TBR) 


iets t * 2 WOonbp IN OUTSIDE BITS IBL & IBR 


it EGER. +2 EN, LO 


LO=IN 


TF (IBR.FQ.0) GO 7) 50 


Pe Site t (iN 18 R) 


fea chee tO, LB l= 15-1 BR} 


LO=IF 


mons KT SEtri(t) ,iSs-1 BL) 


RETU 5N 


END 
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//GO.FT20F001 DD UNIT=3400-4,VIL=SER=GMDT1B, DISP=(OLD, PASS), 


ide LABEL=(1,NL, ,IN), 


77 DCB= (REC FM=FB,LRETL=32, BLKSIZE=512, DEN=2) 


//GO.FT21F001 DD UNIT=3330V,MSV3P=PUBUA, DISP=(OLD,KEEP), 


77 DoN=tS5. 51592.6MDI1A, 


// DC B= (REC FM=VBS,BLKSIZE=4096, LRECL=40 92) 


//GO.-SYSDUME DD SYSOUT=A 


VAGO.SYSIN DD * 


weoem, PSD, MPRY Bay, 1122-1254 LICAL, 17AUG82, 


MmeiGope S25/on0,oV70LT, (DB N&*®*2 VS LOG FREQ) 


im =e FSD, MTRY “BAY,1122-1254 LOCAL, 17AUG82, 


20 AVGS, 32S/SEC,5VOLI, (DB NT**2 VS LOG FREQ) 


STOKES NO. ONE, MTRY BAY, 1122-1254 


TOeAL, V7NUIG82,,. 20 AVS, 32S/SEZ, SVOLT 


SHORES NO, TWO, WIRY BAY, 1122-1254 


MOGAL, 17AUGO2Z, 20 AVSS, 32S/SEL, SVOLT 
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Semes NO. THREE, MTRY BAY, 1122-1254 


been, 17/AUGS2, 20 AVSS, 32S/SEC, SVOLT 


meer > NO. ZERO, MIRY BAY, 1122-1254 LOCAL, 


T7AUG82, 20 AVGS, 32S/SEC, 5VOLI, (LOG VS LOG) 


Pao OF X AND Y, MTRY BAY 1122-12584 


meee, t/AUGS2, 20 AVSS, 32S/SEL, SVOLT 


Samer OF X AND Y COILS, MIRY BAY, 1122-1254 


Poe@aL, 17AUGS82, 20 AVSS, 32S/SEC, SVOLT 


Peer SD, LA MESA,1122-12545LICAL,1/AUG82, 


PemeerGs, S52°5/SEC,10VOLT, (DB NT«*2 VS LOG FREQ) 


Peeormr PSD, LA MESA,1122-1254 LOCAL, 17AU0G82, 


BOM AVNGS, 32S/SEC,10VOLT, (DB NT**2 VS LOG FREQ) 


STOKES NO. ONE, LA MESA, 112201254 


Ben, 17hUG82, 20 AVSS, 32S/SEZ, 10VOLT 


SuekES NO. IWO, LA MESA,. 1122-1254 


BOSAL pe ieAG82,. Z08ATS Syas2S/SEL,. 10VOLT 
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~ 


Swemro NO. THREE, LA HESA, 1122-1254 

Peers, V7AUGE2Z, 20 AVSS, 32S/SECT, I0VOLT 

Seeweo NO. ZERO, LA MESA, 1122-1254 LOCAL, 
WemoctoZ, 20 AVGS, 32S/7SEC, 10VOLT, (LOG VS LOG) 
Pieas OF X AND Y, LA MESA, 1122-1254 

Mmerertb, |7AUGSZ, 20 AVSS, 32S/SE-C, 10VOLT 

eomen OF X AND Y COILS, LA MESA, 1122-1254 
Boern, I/AUG8S2, 20 AVYSS, 32S/SE-Z, 10VOLT 
Pamemee POLARIZATION PSD,MTRY BAY, 1122-1254 
Meee, |17/AUGO2, 20 AVGS, 32S/SEC,5VILT, (DB VS LOG) 
Pimemeac POLARIZATION PSD, LA MESA, 1122-1254 
meen, I/AUGez,20 AVGS, 32S/SEC,1) VOLT, (DB VS LOG) 
Meee CURC POLARIZATION PSD, MIRY BAY, 1122-1254 
Boenb, \G(AUGGZ, 20 AVSS ,32S/SEC,, VOLT, (DB VS LOG) 
Pir CIRC PCLARIZATION PSD, LA WESA, 1122-1254 


meni? AuGle2, 20 AVGS, 32S/SEC, 19VDLT, (DB VS LOG) 





COHER RT CIRC POLARIZATION MTRY BAY/LA MESA, 


ez -i254 LOCAL, t/7AUS82, 20 AVSS, 32S/SEC 


e@leR LEFT CIRC POLARIZATION MIFRY BAY/LA MESA, 


ieee 1254 LOCAL, 1780382, 20 AVSS, 32S/SEC 


/* 


// 
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